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Abstract 

The spin - 3/2 Ising model described by the most general Hamiltonian with an 
up-down symmetry 

-I3H = ]T {JSiSj + KSfs] + LSfs] + y (SiSj + SjSf)} - Sf 

<ij> i 

is investigated on a square lattice. It is shown that this model is reducible to an eight 
- vertex model on a surface in the parameter space spanned by coupling constants 
J, K, L and M. It is shown that this model is equivalent to an exactly solvable free 
fermion model along two lines in the parameter space. Consequently the critical 
behaviour and, in particular, the critical temperature for the second-order phase 
transitions of the model was found exactly. 

PACS number: 05.50+q, 75.10.Hk, 75.40.Cx. 
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The Ising model has been one of the most actively studied systems. In particular, spin 
- 1 and spin - 3/2 Ising models, which present a rich variety of critical and multicritical 
phenomena, are of special interest. 

The spin - 1 Ising model with bilinear (J) and biquadratic (K) nearest- neighbor 
interactions and a single-ion potential (A) is know as the Blume- Emery- Griffiths (BEG) 
model Hfl. The spin - 3/2 Ising model was introduced to explain phase transitions in 
DyVO^ and its phase diagrams were obtained within the mean- field approximation 0. 
Later, this model was used in a study of tricritical properties of a ternary fluid mixture 
The complete phase diagram of the spin - 3/2 Ising model model with L = M = has 
been fully analyzed with the use of two different approaches: mean field and Monte-Carlo 
techniques [||. 

Most studies for two-dimensional lattice mentioned above, are performed for the square 
lattice, very few exact results for spin - 1 and spin - 3/2 Ising models have been carried 
out only for the trivalent lattices, such as the honeycomb lattice. Recently Horiguchi 
|J, Wu and Rosengren and Haggkvist using different theoretical approaches have 
exactly solved the BEG model for the honeycomb lattice in the subspace of interaction 
constants J and K 

exp(K) cosh J = 1. (1) 

The same result was obtained also and for the Bethe lattice ||. 

Wu and Wu || and Kolesik and Samaj [|ID have considered the BEG model in an 



external magnetic field (H) and obtained an exact critical line for all values of H. Recently, 
Lipowski and Suzuki O] and Ananikian and Izmailian [I2"| found the conditions under 



which the spin - 3/2 Ising model on the honeycomb lattice has the same partition function 



as the exactly solvable zero- field spin - 1/2 Ising model. Very recently, Horiguchi |]13 
proposed a general method by which a general spin-S* Ising model in expressed in terms 
of an Ising model of spin ±1 and spin less then S. This exact results provides an excellent 
tool for comparing the accuracy of different approximation schemes mentioned above. 

In the present paper, we have solved exactly the most general spin - 3/2 Ising model for 
the two-dimensional square lattice in the subspace of the four- dimentional space spanned 
by the coupling constants J, K, L and M. We show that this model is reducible to an 
eight-vertex model, for which the exact solution is not known at the present except for 



a few special cases |14|, [15| . Moreover, we have established the equivalence of our model 
on the square lattice with one of the special cases, the so-called "free fermion model", 
which is the eight- vertex model under the "free fermion condition", along two lines in the 
four-dimentional space spanned by the coupling constants J, K, L and M. 

1. We consider the most general spin - 3/2 Ising model with an nearest-neighbor 
interaction and an up-down symmetry, which is described by the following Hamiltonian 

-(3H= USiSj + KSfS] + LS?S* + ^-(S^ + SjSf)} - A £ £/?, (2) 

<ij> i 

where Si = ±|, ±| is the spin variable at site i and < ij > indicates the summation over 
the pairs of nearest-neighbor sites. 
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The equilibrium statistics of the system given by Eq. (0) is determined by the partition 
function Z = J2{ s } ex v(~PH), where (3 = 1/ksT is the inverse temperature and the sum 
is over all spin configurations. 

It is not difficult to show (for details see Ananikian and Izmailian [[Tj|) that if 



{tanh 2 J i = tanh J 2 tanh J 
exp(-4K) = cosh(J 2 - J ) 



(3) 



where 



J = -JH MH L, J 1 = -JH MH L, J 2 = -JH MH L, 

° 4 16 64 ' 4 16 64 ' 4 16 64 ' 

we can write the following identity 

exp {JS X S 2 + KSlS 2 2 + LSfSl + y (SxS 2 3 + S 2 Sf)} = 

= a exp[i?(5 2 + S 2 2 - i)]{l + 4tanh J 5iS 2 exp[ J R (S 1 2 + S 2 2 - i)]}, (4) 

where a — exp(i^/16) cosh J and 

/„r-,x coshJi />t,s tanhJ 2 . . 

exp 2i?) = — : — ^exp - , exp 4flo = — 2 (5) 
cosh J 2 9 tanh J 

Then the partition function defined by Hamiltonian in Eq.(ED can be written as: 



Z=(a e-f) E 1] {l + 4tanhJ ^exp[/2 (S 2 + 5 2 --)]}nexp(-A 5 2 ), (6) 

{s} <ij> 1 i 

where E is the total number of edges, Ao = A — ^R, 7 is the coordination number of 
a lattice, and the first product in Eq.(|6|) is extended over all pairs of neighbouring sites. 
This result is valid for any arbitrary lattice. It should be noted there that Eq.(j3|) is an 
analog of the condition in Eq. ([!]) for the spin - 1 Ising model. 

Thus, we obtain the condition given in Eq.([|), which defines the surface in the space 
spanned by the coupling constants J, K, L and M, where the partition function is written 
in the form of Eq.@. Recently, the equivalence of the spin - 3/2 Ising model on the 
honeycomb lattice with a zero-field spin - 1/2 Ising model on the same lattice in the 



subspace given in Eq. (|3|) have been established [|TT|, [L2| . 

2. In the present section, we investigate the most general spin - 3/2 Ising model on 
the square lattice, described by the Hamiltonian in Eq.(fJ). We show that this model is 
reducible to an eight-vertex model and can be exactly solved on the two nontrivial lines 
in the surface given by Eq.(|3|). 

Now, consider a square lattice (where the coordination number 7 is equal to four) 
composed of N sites (or vertices) and of 2N lattice edges. Then we expand the product 
n<jj> in Eq.(||) and represent graphically each term in the expansion as follows: draw a 
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dotted or solid line over the lattice edge (ij) if the corresponding term in the expansion 
contains the factor 1 or 4 tanh Jo SiSj exp [Ro(Sf + S] - §)]. This leads to eight different 
kinds of configurations, shown in Fig.l, that can occur at a vertex. We assign a Boltzman 
weight {wi} to a vertex having 2n solid and 4 — 2n broken edges. When, the summations 
for all sites are carried out, we obtained the following values for weights to vertices: 



£ t n S?exp[(nR -A )S*] = 

=C 2'- C 2 

where t = tanh J and n is the number of lines with site i as an end-point. 

These facts enable us to rewrite Eq.(^|) in the form where the sum is over all line 
configurations on the square lattice having an even number of lines into each site. For the 
square lattice, this leads to an eight- vertex model shown in Fig.l with the vertex weights 




wi = a, W2 = c and = . . . = wg = b 

where a, b and c are given in Eq.(^). 

Thus, we obtain the exact equivalence 



Z = (a e-^) 2A %,({^}), (9) 
where Z$v ({u>i}) is the partition function of eight-vertex model given by 

Z 8v {{ Wi }) = 2 Yl Ii w i ( 10 ) 



all line configurations % 



The factor of 2 comes from the fact that a reversing of all spins leaves the line config- 
urations unchanged. 

The eight-vertex model on the square lattice plays an important role in the study of 
phase transitions in lattice systems. Unfortunately, except in some special cases jnj |15| 
the behavior of this general model is not known. The exact expression of the free energy 
of this model was first obtained by Fan and Wu |Tj| and by Baxter [|TJ[] in respective 
conditions. In particular, the former authors solved the "free fermion model". This 
model is defined as a particular case of the eight- vertex model in which the vertex weights 
satisfy the relation: 

WiW 2 + w 3 w 4 = w 5 w 6 + w 7 w 8 , (11) 

which is called the "free fermion condition", as in the S - matrix formulation of the 
eight-vertex problem, this condition is equivalent to the consideration of non-interacting 



many-fermion system [16 



It is readily verified, using Eqs.(^) and (|8|), that the "free fermion condition" in Eq. (|TT|) 
is satisfied if, and only if 

exp(4i? ) = \. (12) 
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This equation together with Eqs.(|^) and (|5]) give us the following two nontrivial lines 
in the four- dimensional parameter space spanned by coupling constants J, K, L and M, 
i. e., 

/ (i) 16J = 49L = -14M, K = 0; ( J = £ J) 



(ii) 16J = 169L = -26M, K = 0; ( J = ^ J) 



(13) 



on which the spin - 3/2 Ising model is equivalent to exactly solvable free fermion model. 
Here, we remark that except at the trivial point (J = L = M = K = 0) for special 
(M = 0) spin - 3/2 Ising model, in general the exactly solvable case can not be obtained 

The equivalence, given explicity by Eqs.@, (|10"D and fll3D permits us to deduce exact 
analytic properties of the general spin - 3/2 Ising model. 

The condition given by Eq . (p~3|) can be also interpretated in another way. Consider, 
for example, the case (i): Introducing the two new spin variables <7j and t{ 

*=|s<(S?-£), U = Sf-\ (14) 
we can write the Hamiltonian given by Eq.(||]) in the form 

D <ij> i 

As it easy to see from Eq. (|l4"D the two states Si = |, — \ were transformed into cij = 1 
and two states Si = — |, | into a,; = —1, while the states Si — ±| and Si = ±| into ti — 1 
and tj = — 1 respectively. This gives a one-to-one correspondence between Si and a pair 
of spins (crj,tj). Sowe can say that the our new spin variables a, and £j are independent 
and consequently the summation in J2ia,t) exp(— (3H) can be carry out separatly. 

Now, we consider the thermodynamic properties of our model. A closed expression 



for the free energy of the free fermion model is well known ]TJ| and after some algebraic 
manipulation, we obtain, in the large iV limit, the free energy for the spin - 3/2 Ising 
model on the square lattice in the subspace given by Eq.([L3|) 

-Pf = ln{2v / 3exp(- — )[1 + exp(-2A)]} + — / d& / dy?ln[c 2 + s(cos^ + cos(^)], 

4 87H jo Jo 

(15) 

where c = cosh2J , s = sinh2J and J is determined from Eq.(|l3|). 

Since the second derivative of the free energy in Eq.(^) diverges logarithmically at 
s — s c — 1, the critical behavior of our model are summarized as follows: The spin - 3/2 
Ising model exhibits a first-order phase transition if s > 1, a second-order phase transition 
if s = 1 and no transition at all if s < 1. Thus, a second-order phase transition occurs 
at a temperature determined by sinh2J = 1. This critical condition gives us the two A - 
lines in the space spanned by J and A 

49 r 

(%) J = — ln(l + v 2) = 3, 0886... and A — arbitrary, 
18 
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{%%) J = ln(l + y/2) = 2, 6631... and A — arbitrary, 

where the spin - 3/2 Ising model exhibits an Ising-type phase transition (logarithmic 
specific heat singularity). 

Besides the intrinsic interest surrounding the spin - S Ising model and possible ap- 
plications in real physical situations, one is further attracted to the search for obtaining 
exact nontrivial solutions. So far soluble problems are very few in numbers. In this pa- 
per, an exact solution for most general spin - 3/2 Ising model is obtained. It is shown 
that this model is equivalent to exactly solvable free fermion model along two lines in the 
parameter space given by Eq. (|13"D . In particular, the analytical expressions for the free 
energy per spin the A - lines of Ising-type phase transition were found exactly. 

Finally, we point out, that our result suggests that equivalence of spin - 3/2 Ising 
model to the free fermion model can be extended to a most general Ising model with an 
half-integer spin. That is, the existence of the exact solvable case for the Ising model with 
an half-integer spin results from the absence of the Sf = state. 
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